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ABSTRACT 

Quantitative  information  on  the  response  of  soil  in  natural  deposits 
to  lateral  loadings  was  obtained  by  the  performance  of  field  tests  in  con¬ 
junction  with  theoretical  studies  and  in-situ  and  laboratory  determinations. 
The  field  tests  utilized  a  segmental  pile,  a  lateral  plate  loading  device,  rigid 
poles  subjected  to  lateral  loads,  and  full-scale  laterally  loaded  piles.  Results 
of  the  final  series  of  field  tests  and  analytical  studies  in  the  program  are 
presented  in  this  report,  along  with  the  design  recommendations  formulated 
on  the  basis  of  the  research. 

Procedures  were  developed  for  the  prediction  of  nonlinear  lateral 
soil  pressure-displacement  relationships  applicable  to  undisturbed  deposits 
of  both  cohesive  and  noncohesive  soils  using  easily  measured  conventional 
engineering  properties  of  the  soil.  These  procedures  were  utilized  in  the 
development  of  design  procedures  for  laterally  loaded  piles  in  nonhomo- 
aeneous  layered  soil  systems.  An  iterative  solution  of  the  laterally  loaded 
pile  problem  utilizing  a  digital  computer  was  developed,  and  a  similar 
iterative  solution  for  hand  computations  in  a  tabular  form  was  evolved  for 
use  in  cases  where  a  computer  cannot  be  used. 
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INTRODUCTION 


Subject  and  Purpose  of  Report 

Presented  in  this  report  are  the  results  of  a  study  conducted  at  the 
Naval  Civil  Engineering  Laboratory  (NCEL)  to  develop  realistic  criteria  for 
the  design  and  analysis  of  laterally  loaded,  pile-supported  foundations.  The 
overall  objective  of  the  program  was  to  develop  procedures  for  determining 
soil  moduli  and  variations  of  soil  moduli  with  depth,  width  of  loaded  area, 
magnitude  of  displacement,  and  repetition  of  loading.  For  expediency,  the 
study  was  directed  toward  a  consideration  of  the  laterally  loaded  pile 
problem,  in  order  that  the  relatively  large  body  of  theoretical  and 
experimental  information  on  that  subject  could  be  exploited.  Results  of 
segmental  pile  tests,  lateral  load  tests  on  full-scale  piles,  and  tests  to 
determine  in-situ  soil  properties  were  used  to  fulfill  the  objectives  of  the 
study.  Details  of  these  tests  and  the  associated  results  are  outlined. 

Of  particular  interest  in  the  study  of  laterally  loaded  piles  was  the 
development  of  a  design  method  capable  of  accounting  for  the  nonlinear  soil 
pressure-displacement  characteristics  {p-y  curves)  of  naturally  occurring 
heterogeneous  soil  deposits.  The  procedures  which  were  evolved  to  relate  the 
load-displacement  characteristics  of  soil-pile  systems  utilize  a  closed-form 
expression  comprised  of  parameters  quantifying  the  field  boundary 
conditions.  The  resulting  expression  takes  the  form  of  a  relatively  simple 
rectangular  hyperbola,  and  a  major  portion  of  the  research  effort  has  beers 
devoted  to  a  correlation  of  the  terms  in  the  hyperbola  to  measurable 
physical  properties  of  the  pile  and  the  soil  deposit.  Results  of  the  studies  are 
presented  in  the  form  of  a  digital  computer  solution  or,  as  an  alternate,  a  less 
exact  but  more  expedient  slide  rule-type  of  solution  in  tabular  form. 

Background 

Experimental  Work.  Test  programs  of  several  types  have  been 
performed  at  NCEL  in  order  to  provide  quantitative  information  on  the  p-y 
curves  of  soil.  Rigid  pole  tests,1  iateral  plate  loading  tests,1  • 2 -3  lateral  load 
tests  on  full-scale  piles,4  and  segmental  pile  tests4- 5  have  been  performed  in 
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deposits  of  granular  and  cohesive  soils.  The  data  obtained  have  indicated  the 
form  of  the  resulting  p-y  curves,  and  they  have  given  a  tentative  indication  of 
procedures  for  computing  the  necessary  boundary  conditions  for  the  pile 
solutions.  Subsequent  to  the  time  that  these  tests  were  performed,  additional 
segmental  pile  tests  and  lateral  load  tests  on  full-scale  piles  were  performed 
in  a  soil  deposit  at  El  Cenuo,  California,  in  order  to  provide  further 
background  for  the  prediction  of  p-y  cua'es.  Data  from  these  later  tests  are 
presented  in  this  report. 

Pile  Theory.  The  basic  differential  equation  for  a  vertically 
embedded  pile  subjected  to  a  lateral  load  and/or  a  bending  moment  is 


=  -pB 


(1) 


where  El  is  the  flexural  stiffness  of  the  pile,  y  is  the  horizontal  displacement 
of  a  differential  longitudinal  segment  of  the  pile  located  at  a  depth  Z  below 
the  ground  surface,  the  term  p  represents  the  corresponding  average 
horizontal  soil  pressure  on  the  pile  at  that  depth,  and  B  denotes  the  width  of 
the  pile.  Equation  1  is  applicable  for  cases  where  the  axial  load  on  the  pile  is 
less  than  about  10%  of  the  buckling  load;6  for  greater  axial  loads,  a 
second-order  differential  term  is  required.  The  quantity  pB/y  is  usually- 
referred  to  as  the  modulus  of  horizontal  subgrade  reaction,  k. 

The  complexity  of  Equation  1  has  forced  analysts  to  make 
simplifying  assumptions  regarding  the  interrelationships  between  the  various 
parameters  affecting  pile  response,  and  although  some  of  the  assumptions 
have  been  proven  conclusively  to  be  in  error,  their  use  has  been  continued 
because  of  a  lack  of  a  suitable  substitute.  For  example,  the  nonlinear  p-y 
curves  have  been  represented  by  a  linear  relationship  as  indicated  by  the  dash 
line  in  Figure  1.  Use  of  this  assumption  can  result  in  quite  conservative 
results  for  small  displacements  or  possibly  unconservative  results  for 
displacements  greater  than  y, .  At  any  rate,  it  is  evident  that  the  factor  of 
safety  in  the  analysis  would  vary  considerably  with  the  magnitude  of  the 
displacement.  This  variation  is  especially  prominent  in  the  top  few  feet  of 
soil  where  the  greatest  pile  displacements  occur  and  where  the  most 
resistance  to  horizontal  displacement  is  manifested.7  The  resulting  error 
could  be  significant  for  sea  floor  foundations  or  for  other  installations  in 
sediments  of  low  consistency  where  relatively  large  displacements  might 
occur. 
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y1 

Displacement,  y  (in.) 


Figure  1.  Typical  horizontal  soil  pressure-displacement  diagram  for  a  pile  element. 

Another  simplifying  assumption  which  has  been  used  generally  is 
associated  with  the  variation  in  the  assumed  values  of  k  with  depth  in  the 
soil  deposit.  Solutions  to  the  problem  of  laterally  loaded  piles  have  been 
obtained  with  the  variation  of  k  with  depth  being  assumed  to  be  constant,7 
directly  proportional,8  a  step-function,7  a  hyperbolic  function,9  or  an 
exponential  function.10  These  assumed  relationships  of  k  with  depth  are 
illustrated  in  Figure  2.  It  is  clear,  however,  that  the  deficiencies  associated 
with  these  assumptions  are  a  result  not  only  of  the  typical  nonhomogeniety 
of  naturally  occurring  soil  deposits,  but  also  as  a  result  of  variations  in  the 
magnitudes  of  horizontal  displacement  along  the  length  of  the  pile  and 
nonlinear  p-y  curves.  A  more  rigorous  approach  to  the  problem  is  necessary 
in  order  to  incorporate  the  nonlinear  soil  pressure-deformation-depth 
characteristics  of  naturally  occurring  soil  deposits. 

Recent  advancements  in  computer  technology  have  removed  most  of 
the  restrictions  on  mathematically  solvable  variations  of  soil  .pressure  with 
depth  and  displacement.  Analog  computers  have  been  used  to  study  the 
behavior  of  laterally  loaded  piles,7  and  digital  computers  have  supplied  a 
solution  to  Equation  1  by  means  of  a  difference-equation  method  and  an 


3 


iterative  process  capable  of  adjusting  a  linear  soil  parameter  until  soil-pile 
compatibility  was  obtained.8  Either  type  of  computer  is  capable  of 
providing  a  solution  which  incorporates  nonlinear  soil  characteristics  and 
varying  pile  geometry  for  which  Equation  1  becomes 


f(Z) 


d4y 

dZ4 


-f  (y.z) 


(2) 


where  the  flexural  stiffness  of  the  pile  and  the  nonlinear  soil  pressure- 
displacemept  characteristics  are  both  functions  of  depth.  With  the  aid  of  a 
digital  computer,  it  was  possible  during  this  study  to  consider  all  parameters 
influencing  lateral  soil  pressure-displacement  relationships  without  making 
any  attempt  to  minimize  or  camouflage  the  effects  of  nonlinearity. 


Modulus  of  Horizontal  Subgrade  Reaction,  k  (psi) 


Figure  2.  Assumed  variations  of  k  with  depth. 
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A  graphical  representation  of  the  general  problem  is  presented  in 
Figure  3  for  a  typical  flexible  pile  subjected  to  a  groundline  shear,  Qg,  and  a 
groundline  moment,  Mg.  Under  these  loading  conditions,  the  pile  would 
deflect  approximately  as  shown  in  Figure  3{b).  At  various  locations  along 
the  length  of  the  pile,  the  p-y  curves  might  be  similar  to  those  illustrated  in 
Figure  3(c)  where  the  displacement  corresponding  to  each  location  is 
denoted  by  a  cross.  It  is  apparent  that  any  distribution  of  k  with  depth  could 
result,  and  it  is  equally  evident  that  the  resulting  magnitudes  of  k  will  change 
when  the  pile  loading  is  altered.  Thus,  the  problem  becomes  one  of 
predicting  the  p-y  curves  for  a  given  soil-pile  system. 


(a)  Loading 
conditions. 


(b)  Displacement. 


(c)  Soil  pressure-displacement 
characteristics. 


Figure  3.  Graphical  representation  of  the  problem. 


Rectangular  Hyperbola.  A  technique  for  describing  a  soil's  stress- 
strain  (or  pressure-deformation)  relationship  has  been  proposed.  This 
method,  which  is  based  on  the  characteristics  of  a  rectangular  hyperbola,1 1 
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was  found  to  provide  an  adequate  simulation  of  the  stress-strain  relations 
examined.  Also,  certain  of  the  hyperbolic  parameters  were  found  to  be 
related  to  the  physical  characteristics  of  the  soil.  These  parameters  simply 
describe  a  soil's  nonlinear  stress-strain  relationship  in  linear  form  as  described 
below. 

Consider  a  set  of  coordinate  axes  p'  and  y’  with  a  hyperbola  as  drawn 
in  Figure  4,  with  the  axes  as  asymptotes.  If  the  p'  axis  were  translated  to  the 
rignt  to  the  position  of  the  p  axis  and  if  the  y'  axis  were  translated  down¬ 
ward  to  the  position  of  the  y  axis,  then  the  hyperbola  would  resemble  the 
nonlinear  soil  pressure-displacement  curve  shown  in  Figure  1 .  This  hyperbola 
has  two  properties  that  easily  permit  it  to  be  expressed  by  an  equation  in  the 
linear  form  shown  in  Figure  5.  (Note  that  the  ordinate  of  Figure  5  is  the 
abscissa  divided  by  the  ordinate  of  Figure  4.)  These  two  properties  are:  (1) 
the  inverse  of  the  slope  of  a  tangent  to  the  hyperbola  at  the  translated  origin, 
B/kj,  (Figure  4),  represents  the  intercept  of  the  straight  line  on  the  vertical 
axis  of  Figure  5;  and  (2)  the  inverse  of  the  ultimate  value  of  the  hyperbola  in 
Figure  4,  1/pf,  approximately  represents  the  slope  of  the  straight  line 
resulting  in  Figure  5.  Therefore,  the  equation  of  the  rectangular  hyperbola  is 

_v._J._l  l/_v_\ 

pB  k  k;  +  pf  l  B  /  [  ' 

Thus,  the  nonlinear  p-y  curves  can  be  represented  by  the  linear  equation  of 
the  rectangular  hyperbola  if  the  magnitudes  of  the  parameters  kj  and  pf  are 
known.  Procedures  for  predicting  these  parameters  from  measurable  soil 
properties  and  pile  characteristics  are  presented  later. 

Approach  and  Scope 

Field  tests  with  the  NCEL  12-inch-diameter  segmental  pile  have  been 
performed  in  conjunction  with  lateral  load  tests  on  full-scale  flexible  piles  in 
deposits  of  cohesive  and  noncohesive  soils.  Steel  pipe  piles  with  diameters  of 
approximately  4,  8,  12,  and  16  inches  were  used.  The  results  of  the 
segmental  pile  tests  have  been  used  as  a  guide  in  the  formulation  of  proce¬ 
dures  for  predicting  the  parameters  in  Equation  3.  Results  of  the  lateral  load 
tests  on  piles  have  provided  further  information  on  these  parameters  as  well 
as  a  means  for  checking  the  applicability  of  the  theoretical  work.  In 
addition,  published  data  from  lateral  load  tests  on  piles  have  been  utilized  to 
the  extent  possible  in  the  evaulation  of  the  proposed  design  procedures. 
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Figure  4.  Transposed  coordinate  axis  with  typical  hyperbola. 


Figure  5.  The  rectangular  hyperbola,  a  linear  form  of  the  hyperbola. 
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Results  of  field  tests  m  a  granular  hydraulu  till  and  m  a  deposit  ot 
bav  mud  with  a  desiccated  t  rust  were  reported  earlier  4  however  siqntlu  ant 
information  from  those  earlier  tests  are  referred  to  m  this  report  Results  oi 
a  later  series  of  tests  in  a  deposit  ol  cohesive  soil  near  El  Centro.  California, 
are  presented  also  in  this  report  Data  from  the  segmental  pile  tests,  the 
lateral  load  tests  on  piles,  and  the  m-situ  and  laboratory  tests  for  the 
determination  of  soil  properties  at  the  Ei  Centro  site  have  been  used  in 
conjunction  with  previously  reported  test  results  in  order  to  select  realism 
criteria  for  the  design  of  laterally  loader!  piles  These  criteria  are  outlined  for 
use  in  conjunction  with  a  digital  computer  program  as  listed  in  Appendi  ■  B 
In  addition,  a  simplified  design  method  in  tabular  form  is  proposed  for  \.se 
when  a  computer  is  not  available  This  simplified  procedure  can  be  effected 
with  only  a  slide  rule  and  existing  design  charts  however,  the  same  design 
criteria,  with  nonlinear  p-y  curves,  is  applicable  in  the  simplified  approach  as 
well  as  the  computer  solution  An  example  problem  is  solved  in  order  to 
provide  an  illustration  of  the  proposed  design  procedures 


TEST  PROGRAM  AND  PROCEDURES 
Description  of  Test  Site 

Field  tests  (the  final  series  in  the  program)  were  performed  at  E! 
Centro  Naval  Air  Station,  El  Centro.  California.  That  site  was  chosen  because 
it  provided  a  conveniently  accessible  deposit  of  silty  clay  with  a  desiccated 
zone  near  the  surface.  Tests  were  performed  at  four  prepared  site  locations 
consisting  of  two  locations  with  the  soil  at  its  natural  moisture  content  (dry 
locations)  and  the  other  two  iocations  with  water  standing  at  the  surface 
(flooded  locations).  The  flooding  was  intended  to  reduce  the  effect  of  large 
negative  pore  pressures  in  the  desiccated  zone,  thereby  providing  two 
contrasting  soil  conditions  for  the  investigation.  Flooding  was  accomplished 
by  constructing  a  small  embankment  around  the  area  to  be  tested  and 
maintaining  water  to  a  depth  of  1  inch  for  several  days  before  testing.  No 
tests  were  conducted  until  vane  shear  strength  values  stabilized  in  the 
flooded  soil.  One  of  the  flooded  locations  was  prepared  for  the  segmental 
pile  tests  and  the  other  was  prepared  for  the  lateral  load  tests  on  piles. 

At  each  location,  field  vane  shear  measurements  were  obtained  with 
a  5-inch-long  by  2V2-inch-diameter  vane  which  was  hand-driven  into  the  soil 
with  a  sledge  hammer  and  turned  by  hand  with  a  torque  wrench.  A 
maximum  torque  was  developed  at  approximately  5  seconds,  and  the 
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maximum  torque  value  was  recorded  lor  each  incremental  depth  Torque 
corrections  were  made  for  the  soil  friction  on  the  vane  rod  for  each  vane 
shear  value  measured.  The  resulting  values  are  presented  in  a  later  section  of 
the  report 

Soil  specimens  from  the  upper  10  feet  of  soil  at  each  location  were 
obtained  with  3-inch-diameter,  thin-walled  sampling  tubes  using  a  floating- 
piston  sampler  The  samples  were  returned  to  the  laboratory  where 
classification  tests  and  consolidated-undrained  triaxia!  shear  tests  were 
performed.  Pertinent  results  of  these  tests  also  are  presented  in  a  later 
section  of  this  report,  along  with  the  soil  descriptions. 

Procedures  for  the  Segmental  Pile  Tests 

The  segmental  piie  was  developed  to  provide  a  means  for  studying 
the  lateral  load-supporting  characteristics  of  naturally  occurring,  undisturbed 
deposits  of  soil  This  device  consists  of  12-inch-OD  steel  tubing  installed  in  a 
soil  deposit  in  three  separate  longitudinal  segments  with  all  soil  excavated 
from  the  interior  of  the  segments  The  middle  or  test  segment  can  vary  in 
height  and,  during  the  El  Centro  tests,  segment  heights  of  8  and  12  inches 
were  used  A  test  is  performed  by  forcing  the  middle  segment  horizontally 
while  corresponding  loads  and  displacements  are  recorded.  A  cross  section  of 
the  segmental  pile  prepared  for  installation  with  a  4  inch  test  segment  is 
presented  in  Figure  6a  In  Figure  6b,  the  device  is  shown  after  driving  has 
been  accomplished  In  that  sketch,  the  mandrel  with  its  encased  specimen  of 
soil  has  been  removed,  leaving  the  segmental  pile  with  a  vacated  interior  to 
allow  the  insertion  of  the  loading  mechanism.  Testing  at  a  greater  depth  can 
be  accomplished  by  driving  the  segmental  pile  deeper,  after  realignment  of 
the  segments,  removal  of  the  loading  mechanism,  and  reinsertion  of  the  man¬ 
drel  A  retrieval  system  was  developed  to  recover  the  two  lower  segments. 

With  the  mandrel  in  place,  each  pile  was  driven  into  the  soil  by  a 
drop-weight  type  pile  driver  to  a  depth  such  that  the  center  of  the  test 
segment  was  located  at  least  30  inches  below  the  ground  surface.  At  this 
point,  the  mandrel  and  the  interior  soil  were  removed  and  the  loading 
mechanism  was  inserted  At  the  begmninn  of  each  test,  small  pressure  incre¬ 
ments  of  1  psi  or  less  were  applied  to  the  test  segment  in  order  to  define  the 
initial  portion  of  the  soil  pressure-displacement  curves,  larger  increments 
were  used  throughout  the  remainder  o*  the  tests  Each  pressure  level  was 
maintained  constant  while  displacement  time  measurements  were  recorded 
The  load  level  was  changed  when  the  segment  displacement  became 
negligible  during  a  6-mmute  interval  Only  the  final  soil  pressure- 
displacement  relationship  was  t  onsidered  in  the  analysis  of  the  data. 
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Section  A-A 


Section  B-B 


(a)  During  installation. 


(b)  During  testing. 


Figure  6.  Segmental  pile  and  loading  mechanism. 
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The  magnitude  of  lateral  displacement  was  limited  to  0.5  inch  or  the 
15,000-pound  load  capacity  of  the  loading  mechanism.  On  reaching  the 
capacity  of  the  loading  mechanism,  the  pressure  on  the  test  segment  was 
removed  in  4  or  5  approximately  equal  increments.  The  maximum  previous 
pressure  was  then  reapplied  and  removed  as  many  as  20  times  if  displace¬ 
ment  limitations  permitted.  Displacement  readings  were  taken  at  only  the 
maximum  load  and  after  removal  of  the  load  during  these  repetitions. 
Subsequently,  the  segment  was  realigned,  the  loading  mechanism  was  removed, 
the  mandrel  was  inserted,  and  the  pile  was  driven  to  another  test  depth. 

The  program  for  the  segmental  pile  tests  at  El  Centro  is  provided  in 
Table  1.  After  the  first  three  tests  in  the  dry  area  (D1  through  D3)  had  been 
performed,  it  was  learned  that  misalignment  of  the  loading  mechanism  had 
created  a  binding  of  the  loading  mechanism  with  the  upper  segment, 
resulting  in  unreliable  measurements  during  those  tests.  The  remainder  of  the 
segmental  pile  test  data  were  of  acceptable  quality,  and  these  data  are 
discussed  later. 


Table  1.  Test  Program  for  the  Segmental  Pile  in  El  Centro  Soil 


Test  No. 

Height  of  Test 
Segment 
(in.) 

Depth 

(in.) 

Flooded  Area 

FI 

8 

30 

F2 

8 

60 

F3 

8 

.  108 

F4 

8 

36 

F5 

12 

48 

F6 

12 

78 

F7 

12 

108 

Dry  Area 

D4 

12 

30 

D5 

12 

60 

D6 

12 

90 

D7 

12 

114 

D8 

8 

36 

D9 

8 

72 

DIO 

8 

110 

Dll 

8 

36 

D12 

12 

42 
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Procedures  for  the  Lateral  Load  Tests  on  Piles 


Lateral  load  tests  on  flexible  piles  were  conducted  to  provide 
information  on  pile  response  for  use  in  a  theoretical  study.  A  range  of  pile 
sizes  was  considered  in  order  to  enable  an  evaluation  of  scaling  effects. 
Open-ended,  steel,  pipe  piles  of  four  different  sizes  ranging  from  4.5  to 
16  inches  in  diameter  were  driven  to  depths  such  that  a  free-standing  portion 
of  approximately  3  feet  remained  for  each  pile.  The  piles  were  designed  to 
have  embedded  lengths  greater  than  4R  where  R  is  a  relative  stiffness  factor 
defined  as 


R 


dir 


(4) 


The  quantity  kQ  isavalueof  the  modulus  of  horizontal  subgrade  reaction 
for  a  case  where  the  modulus  is  assumed  to  be  constant  with  depth.  It  has 
been  shown  analytically  that  piles  with  embedded  lengths  greater  than  about 
4R  respond  to  lateral  loadings  approximately  as  if  they  were  embedded  to  an 
infinite  depth.7  The  impropriety  of  assuming  that  k  can  be  expressed  as  a 
simple  function  of  depth  without  considering  magnitudes  of  displacement 
and  other  boundary  conditions  has  been  discussed  previously,  but  the  use  of 
the  assumption  for  computing  necessary  pile  lengths  is  an  acceptable 
procedure  if  conservatively  low  values  of  k0  are  assumed.  A  photograph  of 
the  piles  embedded  at  a  flooded  location  is  presented  in  Figure  7. 

During  the  tests,  lateral  loads  were  applied  to  each  pile  at  a  point 
approximately  32  inches  above  the  ground  surface  in  order  to  provide  both  a 
shearing  force  and  a  bending  moment  in  the  pile  at  the  groundline.  A  typical 
test  setup  is  shown  in  Figure  8  where  the  cable  and  hydraulic  ram  used  for 
applying  the  loads  can  be  observed.  Also  shown  (less  distinctly)  in  Figure  8 
are  the  four  dial  indicators,  spaced  approximately  evenly  up  to  30  inches 
above  the  ground  surface,  for  measurements  of  displacement  and  slope,  and 
the  Dillon  dynamometer  used  to  obtain  load  measurements. 

The  program  for  the  lateral  load  tests  on  piles  is  given  in  Table  2 
where  the  test  numbers,  pile  diameters,  pile  stiffnesses,  and  depths  of 
embedment  are  presented.  Each  lateral  load  test  was  conducted  in  a  manner 
similar  to  the  segmental  pile  tests.  The  lateral  loads  were  applied  incre¬ 
mentally  up  to  a  maximum  of  20  kips  (the  capacity  of  the  loading  system) 
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or  until  the  range  of  the  two  top  dial  indicators  was  exceeded  (about  1.5 
inches  of  displacement  at  the  groundline).  Each  load  was  again  maintained 
constant  until  progressive  displacement  during  a  5-minute  interval  became 
negligible. 


Table  2.  Program  for  Lateral  Load  Tests  on  Piles  in  El  Centro  Soil 


Test 

Pile 

Pile 

Depth  of 

No. 

Diameter 

Stiffness 

Embedment 

(in.) 

(lb-in.2) 

(in.) 

Flooded  Area 

PI 

4.50 

2.26  x  108 

144 

P2 

8.625 

1.73  x  109 

204 

P3 

12.75 

8.38  x  109 

264 

P4 

16.00 

1.686  x  10’ 0 

324 

Dry  Area 

P5 

4.50 

2.26  x  108 

144 

P6 

8.625 

1.73  x  109 

204 

P7 

12.75 

8.38  x  109 

264 

P8 

16.00 

1.686  x  10* 0 

324 

3 


Figure  8.  Lateral  load  test  in  progress. 
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TEST  RESULTS 


Data  obtained  from  the  soil  property  determinations,  the  segmental 
pile  tests,  and  the  lateral  load  tests  on  piles  is  presented  and  discussed.  An 
analysis  of  the  test  data  is  included  in  a  iater  section  of  the  report. 

Soil  Properties 

The  soil  at  the  El  Centro  site  consists  of  a  brown,  silty  clay  and 
clayey  silt  with  an  occasional  trace  of  fine  sand.  It  is  classified  as  ML  to  CL 
by  the  unified  soil  classification  system.  A  desiccated  crust  approximately  60 
inches  in  depth  exists  at  the  site,  and  the  water  table  was  located  at  a  depth 
of  approximately  144  inches. 

Index  properties  of  the  soil  are  exhibited  in  Figure  9,  and  the  field 
vane  shear  strength  test  results  are  plotted  in  Figure  10.  Values  plotted 
represent  averaged  values  for  several  tests  at  each  location,  and  the  test  data 
for  the  shallow  depths  of  interest  typically  fell  within  plus  or  minus  30%  of 
the  values  plotted.  It  is  noted  that  the  effect  of  flooding  apparently  becomes 
negligible  at  depths  below  50  inches.  In  the  zone  from  120  to  160  inches 
depth,  the  fairly  large  difference  in  the  values  obtained  for  the  two  areas  is 
believed  to  be  a  result  of  the  lack  of  resolution  of  the  measurements, 
brought  about  by  the  relatively  large  corrections  for  rod  friction  which  were 
necessary  at  the  greater  depths.  The  discrepancy  is  of  no  consequence 
because  the  behavior  of  the  laterally  loaded  piles  is  controlled  almost 
entirely  by  the  characteristics  of  the  soil  at  shallower  depths. 

The  effect  of  desiccation  on  the  vane  shear  strengths  is  readily 
apparent  in  Figure  10.  It  can  be  seen  that  the  reductions  in  negative  pore 
pressures  in  the  desiccated  ctust  as  a  result  of  flooding  led  to  appreciably 
lower  shear  strengths  in  the  flooded  zone.  Therefore,  considerably  greater 
pile  displacements  for  corresponding  magnitudes  of  loading  were  expected 
for  the  flooded  location  in  comparison  with  the  dry  location. 

Segmental  Pile  Test  Data 

When  soil  pressure  was  plotted  versus  corresponding  magnitudes  of 
displacement,  most  of  the  data  resulted  in  S-shaped  curves,  i.e.,  a  dispropor¬ 
tionately  large  amount  of  displacement  occurred  at  low  pressure  levels.  The 
data  points  from  one  typical  test  are  presented  in  Figure  1 1 .  The  reversed 
curvature  resulted  from  alignment  difficulties  encountered  with  the  loading 
mechanism  at  the  start  of  the  tests.  Consequently,  a  method  for  locating  the 
correct  zero  displacement  was  needed,  and  two  methods  were  employed. 
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The  first  method  involved  simply  correcting  the  initial  portion  of  the  curve 
by  eye.  A  second  method  involved  the  plotting  of  soil  pressures  on  an 
arithmetic  scale  versus  corresponding  magnitudes  of  displacement  on  a 
logarithmic  scale.  The  upper  portion  of  the  S-shaped  curves,  which  were 
assumed  to  be  the  most  valid  portion,  plotted  as  a  straight  line;  so  by 
extending  this  straight  line  to  the  logarithmic  axis,  it  was  possible  to  locate 
an  apparent  point  for  zero  displacement.  Both  methods  produced  similar 
results  and  an  average  was  used  when  a  discrepancy  occurred. 


Soil 
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Figure  9.  Index  properties  of  El  Centro  soil. 
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Figure  11.  Typical  S-shaped  soil  pressure-displacement  curve  from 
segmental  pile  tests  with  correction. 

Figure  12  contains  the  corrected  data  obtained  from  segmental  pile 
tests.  The  data  from  the  first  loading  during  each  test  are  presented  in  the 
form  of  soil  pressure-displacement  curves  grouped  according  to  segment  size 
and  test  location.  Each  curve  is  characteristically  concave-downward  and 
similar  to  the  modified  hyperbola  presented  in  Figure  3. 

Data  from  subsequent  cycles  of  the  maximum  load  during  the  tests 
are  presented  in  Figure  13  where  displacement,  expressed  as  a  percentage 
increase  over  the  peak  displacement  occurring  during  the  first  loading  cycle, 
is  plotted  as  a  function  of  the  number  of  load  applications.  An  extrapolation 
of  some  of  the  test  data  to  an  infinite  number  of  repetitions  indicates  that 
more  than  a  100%  increase  in  displacement  as  a  result  of  repetitive  lateral 
loading  in  the  El  Centro  soil  can  be  expected  in  some  cases.  However,  during 
a  majority  of  the  tests  it  appeared  that  the  increase  would  not  have  exceeded 
50%. 
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Figure  13,  Effect  of  repetitive  loadings  on  displacements. 
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It  is  noted  that  the  segmental  pile  test  data  presented  in  Figure  12 
are  not  presented  in  a  form  suitable  for  application  to  the  general  problem  of 
laterally  loaded  piles  It  is  necessary  to  apply  scaling  relationships  and  shape 
(actors  which  have  been  explained  previously  4  These  adjustments  are  dis¬ 
cussed  m  a  later  section  of  the  report  where  the  segmental  pile  test  data  are 
used  as  a  guide  in  the  formulation  of  generally  applicable  nonlinear 
horizontal  load-displacement  relationships  for  soil. 

Data  from  Lateral  Load  Tests  of  Piles 

Lateral  load-displacement  curves  for  the  piles  tested  at  El  Centro  are 
presented  in  Figure  14,  and  similar  curves  for  the  slopes  at  the  groundline  are 
presented  in  Figure  15  It  is  again  pointed  out  that  the  lateral  load  was 
applied  at  a  point  approximately  32  inches  above  the  groundhne  so  that  the 
loading  could  be  resolved  to  an  equivalent  groundline  moment  in  addition  to 
a  groundline  shearing  force  The  effects  of  both  pile  size  and  soil  condition 
can  readily  be  seen  in  the  plots,  and  these  effects  along  with  the  theoretical 
results  illustrated  by  the  solid  lines  in  Figures  14  and  15,  are  discussed 
further 


THEORETICAL  STUDY 

The  theoretical  portion  of  this  program  has  been  associated  mainly 
with  the  development  of  procedures  for  determining  the  nonlinear  lateral 
load-displacement  relationships  for  natural  soil  deposits.  It  has  been  found 
that  a  rectangular  hyperbola  such  as  that  described  earlier  provides  a  satisfac¬ 
tory  representation  of  these  relationships.4  Therefore,  the  following 
discussion  is  devoted  to  the  development  of  techniques  for  relating  the  soil 
properties  and  other  boundary  conditions  to  the  parameters  necessary  for 
the  definition  of  a  hyperbolic  representation  of  the  lateral  load-displacement 
curves 

For  the  definition  of  these  curves  there  are  two  extremely  important 
reference  points  which  must  be  considered.  The  first  of  these  is  the  initial 
slope,  kj,  of  the  curves.  Guidance  in  the  determination  of  that  quantity  can 
be  obtained  from  constant  values  of  the  modulus  of  horizontal  subgrade 
reaction,  kc,  for  preloaded  cohesive  soils  which  have  been  proposed  for  use 
in  predicting  moments  in  laterally  loaded  piles.1 2  The  values  proposed  are 
based  upon  the  consistency  of  the  soil  and  are  approximately  equal  to  the 
product  of  67  times  the  undrained  shear  strength  of  the  soil.  Because  the 
proposed  values  are,  by  necessity,  somewhat  conservative,  and  because  the 
load-displacement  curves  characteristically  exhibit  concave  downward  curva¬ 
ture,  it  is  certain  that  the  initial  moduli,  kj,  will  be  larger  than  67S. 


21 


Figure  14.  Displacements  at  the  groundline  during  pile  tests  at  El  Centro. 
(Solid  lines  represent  computer  solutions.) 


o 


in 


jejajei 


The  other  reference  point  c;in  be  obtained  by  contouring  me  lateral 
bearing  capacity  of  the  soil,  p(,  to  the  vertical  bearing  capacity  at  corres¬ 
ponding  depths.  The  lateral  bearing  capacity  would  be  expected  to  be  no 
larger  than  the  vertical  bearing  capacity  at  large  depths,  and  normally  p, 
would  be  considerably  smaller  than  the  corresponding  vertical  values  at 
locations  near  the  ground  surface  because  ot  the  orientation  of  the  respective 
failure  surfaces.  Thus,  the  ratio  of  p(  to  the  undrained  shear  strength  of  the 
cohesive  soil  would  be  expected  to  be  a  function  of  depth  and,  from  conven¬ 
tional  foundation  engineering  practice,  the  ratio  would  not  exceed  9.0  at 
large  depths  or  5.5  at  shallow  depths. 

With  these  two  reference  points  in  mind,  the  El  Centro  test  data  were 
analyzed  to  determine  the  proper  magnitudes  of  k;  and  p(  for  the  soil  at  that 
location.  Prior  to  analyzing  the  data,  though,  it  was  necessary  to  apply  the 
appropriate  factors  to  account  for  the  size  and  the  shape  of  the  test  segments, 
in  order  to  have  the  segmental  pile  test  data  correspond  to  the  boundary 
conditions  in  existence  during  the  lateral  load  tests  on  piles.  These  factors 
have  been  discussed  in  detail  in  an  earlier  report.4  Summarizing  the  procedure, 
it  is  appropriate  to  plot  1.5d/pHaasa  function  of  1.5d/Ha  where  H  repre¬ 
sents  the  height  of  the  rest  segment,  p  is  the  change  in  the  average  horizontal 
soil  pressure  on  the  test  segment  as  a  result  of  the  corresponding  lateral 
displacement,  d,  and  a  is  a  shape  factor  defined  as  a  =  18  in./(12  in.  +  0.5H) 
for  cases  where  H,  expressed  in  inches,  is  equal  to  or  less  than  12.  It  is  noted 
that  the  quantity  1.5  in  the  above  expression  is  the  assumed  magnitude  of  a 
for  a  laterally  loaded  pile.  With  the  data  expressed  in  this  form,  it  is  possible 
to  determine  the  appropriate  values  of  kj  and  p(  from  the  ordinate  intercepts 
and  the  slopes  of  the  plots,  respectively. 

Initial  Slope  of  Soil  Pressure-Displacement  Curves 

Hyperbolic  plots  of  the  segmental  pile  test  data  were  made  and,  in 
general,  the  plots  were  linear  so  that  magnitudes  of  kj  and  pf  could  be 
calculated.  There  was  considerable  scatter  in  the  values  of  kj  obtained, 
however,  it  was  obvious  that  the  magnitudes  of  kf  were  considerably  larger 
for  the  El  Centro  soil  than  for  the  bay  mud  tested  in  an  earlier  program.4  In 
addition,  there  were  large  differences  in  the  magnitudes  of  k,  between  the 
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flooded  end  the  dry  locations  at  El  Centro,  and  the  magnitudes  of  the 
differences  could  not  be  explained  solely  on  the  basis  of  differences  m 
shearing  strength  of  the  soil  at  the  two  locations  Therefore,  it  was  necessary 
to  do  a  considerable  amount  of  searching  to  determine  the  causes  for  the 
differences. 

In  this  search,  it  was  necessary  to  return  to  the  data  from  the  tests  in 
the  bay  mud.4  Both  the  segmental  pile  tests  and  the  lateral  load  tests  on 
piles  in  that  soil  were  reanalyzed.  It  was  found  that  the  magnitudes  of  k,  for 
that  soil  were  on  the  order  of  100  times  the  corresponding  values  of 
undrained  shear  strength.  The  factor  of  100  is  logical  in  view  of  the  lower 
limit  of  67S  discussed  earlier. 

Returning  to  the  El  Centro  data,  a  similar  study  of  the  hyperbolic 
plots  revealed  that  k,  was  approximately  equal  to  200S  for  the  flooded  area 
and  400S  for  the  dry  area.  It  is  pointed  out  that  the  determination  of  these 
factors  was  not  a  straightforward  matter  of  picking  values  from  the 
hyperbolic  plots.  There  was  such  a  large  amount  of  scatter  in  the  values  of 
k;/S,  resulting  from  the  difficulty  of  experimentally  determining  the  initial 
slope  of  a  lateral  load-displacement  curve  in  soil,  that  it  was  required  to 
consider  k,  simultaneously  with  the  magnitudes  of  p(.  It  should  be  noted 
that  the  measured  values  of  pf  are  relatively  precise  because  they  are 
independent  of  displacement,  whereas  the  measurements  of  k,  require 
accuracy  in  both  soil  pressure  and  displacement.  By  working  backward  from 
known  values  of  p(,  it  was  possible  to  select  values  of  ks/S  which  were 
compatible  with  the  measured  test  data.  In  that  manner,  the  values  of  200 
and  400  were  determined,  and  these  values  were  found  to  be  within  the 
range,  but  conservatively  below  the  mean,  of  the  values  of  k,  computed  from 
the  hyperbolic  plots. 

It  was  found  during  the  study  of  k,  that  the  parameters  of  the 
rectangular  hyperbola,  k,  and  pf,  could  not  be  determined  precisely  from  the 
test  data.  However,  u  is  pertinent  that  in  cases  where  k,  was  found  from  a 
plot  to  be  disproportionately  laige,  the  corresponding  value  of  pf  was  found 
to  be  disproportionately  low.  As  a  result,  the  use  of  the  faclois  200S  and 
400S  for  determining  k,,  along  with  the  use  of  the  corresponding  expression 
for  p(  (presented  later),  provided  a  realistic  representation  of  the  measured 
soil  pressure-displacement  curves. 

It  is  pointed  out  that  computed  magnitudes  of  pile  response  are  not 
very  sensitive  to  changes  in  the  assumed  magnitudes  of  k|(  therefore,  it  is  not 
important  that  k,/'S  be  determined  accurately.  The  reason  for  that  lack  of 
sensitivity  is  suggested  in  Figure  16  where,  for  a  constant  value  of  p(,  the 
limited  effect  of  a  relatively  large  peicentage  change  of  k,/S  on  the  resulting 
soil  pressure-displacement  curves  can  be  seen.  A  demonstration  of  the 
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relatively  small  influence  of  changes  in  k;/S  appears  in  Figure  1 7  where 
theoretical  predictions  of  pile  displacements  are  presented  for  two  values  of 
the  ratio.  It  can  be  observed  that  a  doubling  of  kj/S  resulted  in  only  about  a 
22%  decrease  in  the  displacements  computed. 

Possible  reasons  for  the  resulting  differences  in  the  magnitudes  of 
kj/S  were  investigated.  It  was  suspected  that  the  ratio  was  related  in  some 
manner  to  the  presence  of  negative  pore  pressures  in  the  desiccated  crust,  so 
a  readily  determinable  measure  or  index  of  these  pore  pressures  was  sought. 
It  v/as  feit  that  the  liquidity  index  should  exhibit  a  gross,  qualitative-type  of 
relationship  to  the  negative  pore  pressures,  so  a  plot  of  liquidity  index,  Lj 
versus  k;/S,  Figure  18,  was  made.  Since  neither  kj/S  nor  Lj  could  be  defined 
with  precision,  the  data  were  represented  as  shaded  zones  rather  than  data 
points.  An  apparent  discrepancy  exists  in  the  data  for  the  bay  mud  in  that  a 
change  in  Lf  apparently  caused  no  change  in  k;/S.  No  attempt  is  made  to 
explain  this  discrepancy,  however,  in  view  of  the  general  recommendations 
to  follow,  the  discrepancy  is  inconsequential. 


Segment  Displacement,  d  (in.) 


Figure  16.  Soil  pressure-displacement  curves  for  different  soils  with 
similar  failure  pressures  and  different  initial  slopes. 


Figure  17.  Typical  force-pile  displacement  curves,  emphasizing  the  effect 
of  varying  kj/S  in  computer  solution. 

Based  upon  the  results  of  the  segmental  pile  tests  and  the  lateral  load 
tests  on  piles,  it  is  recommended  that  the  magnitude  of  kj/S  be  taken  as  100 
for  cohesive  soils 


k,  *  100S 


(5) 


27 


It  is  recognized  that  for  desiccated  clays  with  a  low  liquidity  index  the 
recommendation  will  result  in  overconservative  values  of  computed 
displacements;  however,  the  test  results  have  indicated  that  the  value  of  100 
would  be  approached  in  such  a  case  if  the  soil  were  to  become  saturated. 
Thus,  the  use  of  a  value  larger  than  100  could  prove  to  be  unconservative  if  a 
heavy  rainstorm  or  some  other  phenomenon  causing  flooding  should  occur. 
The  use  of  a  value  lower  than  100  is  not  warranted  because  available  test 
data  indicate  that  no  lesser  values  would  occur,  and  that  value  is  realistic  in 
comparison  with  the  lower  limit  of  67  discussed  previously. 


Stiffness  -  Vane  Strength  Ratio,  K ■  /S 


Figure  18.  Relationship  between  hyperbolic  parameter  and  liquidity  index. 
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Lateral  Bearing  Capacity 


During  the  determination  of  an  appropriate  relationship  for  the 
lateral  bearing  capacity,  pf ,  of  cohesive  soil,  it  was  again  necessary  to  return 
to  the  data  pertaining  to  the  bay  mud,  discussed  previously  in  Reference  4. 
Results  of  theoretical  studies  of  lateral  bearing  capacity13  were  used  as  a 
guide  in  the  formulation  of  the  relationship  for  p{  proposed  in  Reference  4. 

It  was  found  that  p{  was  directly  proportional  to  S  and  Z/B  where  S  is  the 
undrained  shear  strength,  Z  is  the  depth  below  the  ground  surface,  and  B  is 
the  pile  width.  An  investigation  of  the  El  Centro  data  indicated  that  the 
same  expression  for  pf  provided  predictions  which  were  within  about  ±30% 
of  the  measured  values.  However,  it  was  later  found  to  be  necessary  to  alter 
the  scaling  relationship,  Z/B,  on  the  basis  of  the  data  from  the  lateral  load 
tests  on  piles. 

The  relationship  for  p{  proposed  in  Reference  4  led  to  realistic 
predictions  of  pile  response  in  the  bay  mud;  however,  it  was  noted  that  a 
straight  Z/B  scaling  resulted  in  values  of  pf  which  were  too  low  for  16-inch- 
diameter  piles  and  slightly  high  for  4.5-inch-diameter  piles.  A  similar 
situation  was  discovered  with  the  El  Centro  data.  Therefore,  the  theoretical 
studies  were  consulted,  and  it  was  noted  that  the  assumption  of  inverse 
proportionality  of  pf  to  B  resulted  in  the  pile  width  having  more  influence 
than  was  appropriate.  Consequently,  other  functions  of  pile  width  were 
investigated,  and  it  was  found  that  the  square  root  of  B  was  appropriate  for 
the  data  from  the  lateral  load  tests  in  both  the  bay  mud  and  the  El  Centro 
soil.  Thus,  the  following  relationship  for  pf  was  evolved 

Pi  *  0.25  S^j  (6) 

It  is  noted  that  Equation  6  is  not  dimensionally  correct.  The  ramifications  of 
the  nondimensionality  have  been  studied,  and  the  conclusion  was  reached 
that  the  resulting  equation  is  only  an  empirical  approximation  of  a  perhaps 
more  complicated  dimensionally  correct  relationship  for  pf.  It  does, 
however,  provide  an  expedient  means  of  computing  fairly  reliable  values  of 
lateral  bearing  capacity  for  cohesive  soils. 

One  further  point  which  must  be  mentioned  in  connection  with  the 
expression  for  p{  is  the  reference  point  alluded  to  earlier  in  regard  to  the 
maximum  possible  value  of  pf/S.  From  Equation  6  it  is  obvious  that  the 
stipulation  that  pf/S  at  shallow  depths  could  not  exceed  approximately  5.5 
is  satisfied.  As  one  would  intuitively  expect  from  the  orientation  of  the 
failure  surfaces,  that  ratio  would  be  far  less  than  the  corresponding  value  for 
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vertical  bearing  capacity.  Also,  it  was  mentioned  earlier  that  pf/S  could  not 
exceed  approximately  9.0  at  great  depths.  This  latter  criterion  is  not  met  by 
Equation  6,  however,  the  values  of  Z/^/B  for  which  p,/S  exceeds  9.0  are  in  a 
zone  of  soil  sufficiently  deep  that  it  has  only  a  negligible  influence  on  the 
response  of  laterally  loaded  piles.  Therefore,  there  was  no  justification  in  this 
study  for  a  more  elaborate  relationship  or  perhaps  a  step  function  for  pf. 

The  criteria  presented  previously  have  been  applied  to  a  large  number 
of  cases  reported  in  the  literature.  In  all  cases,  the  proposed  criteria  have  led 
to  conservative  results  and,  in  some  cases,  the  computed  response  proved  to 
be  somewhat  overconservative.  However,  there  was  no  indication  that  the 
criteria  should  be  altered  to  account  for  these  conservative  results  because 
the  soil  conditions  repoited  in  the  literature  were  not  defined  as  accurately 
as  would  have  been  desirable,  especially  at  the  very  shallow  depths. 
Frequently,  the  soil  data  were  in  the  form  of  qualitative  word  descriptions  of 
the  consistency  of  clay,  e.g.,  soft,  very  soft,  medium,  etc,  or  similar  descrip¬ 
tions  of  the  relative  density  of  sand,  e.g.,  loose,  medium,  etc.  These 
descriptions  were  extrapolated  to  assumed  quantitative  values  of  shear 
strength  parameters,  and  the  resulting  errors  were  apparently  large. 

Other  factors  which  contributed  to  the  conservative  nature  of  the 
computed  results  are  associated  with  the  methods  used  to  determine  the 
reported  soil  properties  and  with  the  tune  at  which  the  properties  were 
measured.  Many  of  the  values  of  consistency  of  the  cohesive  soil  were  based 
upon  the  results  of  unconfmed  compression  tests  performed  with  soil 
samples  obtained  with  thick-walled  samplers.  Therefore,  the  soil  samples 
were  undoubtedly  disturbed.  Other  values  of  consistency  of  cohesive  soils 
were  based  upon  the  results  of  standard  penetration  tests,  an  unreliable 
procedure  for  that  type  of  soil. 

The  effect  of  the  time  at  which  the  soil  properties  were  determined  is 
relevant  only  for  the  tests  in  granular  soil.  It  was  decided  during  the  experi¬ 
mental  program  that  the  prudent  approach  would  be  to  measure  the 
properties  of  the  granular  soil  subsequent  to  the  pile  driving  operations  at 
the  site  because  the  pile  driving  would  tend  to  densify  the  soil,  leaairg  to  a 
stiffer  response  of  the  pile.  Therefore,  if  correlations  were  made  on  the  basis 
of  soil  properties  determined  prior  to  driving,  an  unconservative  design 
procedure  could  result.  On  the  other  hand,  the  approach  taken  may  lead  to 
conservative  results,  if  the  soil  properties  are  measured  prior  to  pile  driving. 
For  the  tests  reported  in  the  literature,  the  properties  were  measured  prior  to 
driving  so  conservative  resuits  were  computed. 
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Computation  of  Lateral  Soil  Pressure-Displacement  Diagrams 

Earlier  discussions  have  described  the  procedure*  used  to  determine 
the  parameters  required  to  define  lateral  soil  pressure-displacen  unt  mi.mon 
ships  m  terms  of  soil  and  pile  properties.  These  relationships  are  summarized 
in  the  following  discussion,  and  similar  relationships  for  granular  soil 
deposits  are  presented. 

Cohesive  Soil.  For  a  determination  of  the  lateral  soil  pressure- 
displacement  relationships  applicable  to  cohesive  soils,  it  is  neressary  to 
obtain  a  measurement  of  the  undrained  shearing  strength  to  i  depth  of 
perhaps  20  pile  diameters.  These  strength  values  can  then  be  substituted  m 
Equations  5  and  6,  repeated  below,  to  determine  the  necessary  parameters 
for  computing  the  relationships 


k,  *=  100S 


(5) 


Pf 


(6) 


These  parameters  can  in  turn  be  substituted  into  Equation  3  to  provide  a 
convenient  representation  of  the  soil  pressure-displacement  relationship 


1  =  X  =  1  +  _l/v\ 

k  pB  k,  pf  \B  / 


(3) 


Granular  Soil.  For  granular  soil,  the  computational  procedures  are 
somewhat  more  complicated  than  the  corresponding  procedures  for  cohesive 
so1  Is.  First,  it  is  necessary  to  obtain  a  measure  of  the  angle  of  internal 
friction,  0,  of  the  soil  to  a  depth  of,  again,  perhaps  20  pile  diameters.  This 
can  be  obtained  by  triaxia!  tests  or  by  index  measurements  of  relative 
density.  During  earlier  tests  in  a  hydraulic  fill  of  granular  material,  it  was 
found  that  the  vane  resistance,  S,  in  psi  (computed  as  if  it  were  vane  shear 
strength)  was  approximately  equal  to  the  standard  penetration  resistance,  N, 
in  blows  per  fool  4  Thus,  the  following  relationship  for  0  was  proposed  on 
the  basis  of  available  data  on  standard  penetration  tests  14- 15 


l  25.4  -  p0  \ 
‘  1 2.94  +  0.1225  p0j 


S0' 


lib 


(S*  N) 


(7) 
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1 1  is  noted  thnt  <p  is  in  decrees  when  S  is  in  psi  mid  when  pQ,  the  olfecti  /o 
overburden  pressure,  is  in  psi. 

In  the  earlier  work,4  the  relationship  ol  4>  to  p,  was  based  entirely 
upon  theory  because  the  segmental  pile  tests  could  not  be  carried  to  failure 
of  the  soil  in  bearing.  However,  subsequent  comparisons  of  computed  and 
measured  pile  response  have  indicated  that  the  proposed  relationships  were 
of  approximately  the  correct  form,  but  the  predicted  magnitudes  of  p,  at 
shallow  depths  were  too  small.  Therefore,  the  original  form  of  the  relation¬ 
ship  was  retained,  and  the  coefficients  were  altered  empirically  in  order  to 
provide  more  realistic  predictions  of  lateral  bearing  capacity  at  shallow 
depths.  The  following  relationship  resulted 


Pf 


Po 


0.664  (IQ0038  <t>)  +  0.04  9  (10  0  034  0) 


(8) 


Here,  as  before,  0  is  expressed  in  degrees  and  p0  is  the  effective  overburden 
stress. 

It  was  found  during  the  previously  reported  work  that 


k;  «  2,000  psi  (9) 

for  sand.4  An  analysis  of  data  from  lateral  load  tests  on  piles  reported  in  the 
literature  has  indicated  that  Equation  9  provides  reasonable  results.  There¬ 
fore,  the  lateral  soil  pressure-displacement  relationships  for  granular  soils  can 
be  computed  with  Equation  3  and  the  input  parameters  provided  by 
Equations  7.  8,  and  9. 

Repetitive  Loadings.  The  effects  of  repetitive  and  cyclic  loadings  can 
be  accounted  for  by  adjusting  the  input  soil  parameters  or  by  applying 
correction  factors  to  the  computed  pile  response.  The  latter  procedure  has 
been  chosen  because  the  necessary  adjustments  to  the  soil  parameters  could 
lead  to  erroneous  predictions  of  moments  in  '.he  piles. 

Cyclic  or  repetitive  loadings  have  the  effect  of  increasing  a  pile 
displacement  at  a  decreasing  rate  with  increasing  numbers  of  cycles  or 
repetitions.  This  results  from  either  remolding  an  undisturbed  cohesive  soil 
or  densifying  a  granular  soil  by  the  process  of  repeated  shearing 
deformations,  The  rosults  of  field  and  model  tests,  os  summarized  by 
Prakash,16  indicate  a  20  to  50%  increase  in  groundline  displacement  for 
infinite  repetitions  of  a  constant  load  in  granular  and  cohesive  soils.  The 
cohesive  soils  tested,  however,  were  insensitive  clays.  Little  is  known  of  the 
effects  of  repetitive  pile  loads  on  sensitive  clays  and  a  discussion  thereof  is 
excluded.  More  research  is  needed  in  the  area. 
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Repetitive  load  tests  were  performed  with  the  segmental  pile,  and  the 
results  of  these  tests  were  shown  m  Figure  13.  Similar  results  from  earlier 
tests  with  the  segmental  pile  in  bay  mud  have  been  reported.4  These  latter 
data  indicate  percentage  increases  in  the  vicinity  of  the  upper  limit  of  the 
range  proposed  by  Prakash,  so  it  is  recommended  that  a  50%  increase  of 
displacements  be  assumed  for  conditions  of  cyclic  or  repetitive  loadings. 

An  investigation  of  the  effects  ol  repetitive  and  cyclic  loadings  on 
the  magnitudes  of  maximum  moments  in  the  piles  indicated  that  there  was 
negligible  danger  of  a  moment  increase  in  a  pile  during  these  additional 
loadings. 

Computed  Pile  Response 

Computational  Procedure  Equation  2  was  solved  by  an  iterative 
finite  difference  procedure  explained  in  an  earlier  report.4  For  this  study, 
the  computer  program  was  refined  somewhat,  and  the  lateral  soil  pressure- 
displacement  parameters  presented  in  the  preceding  sections  were 
incorporated,  1  he  proposed  final  version  of  the  computer  program  is 
presented  in  the  Appendix 

For  each  computation,  the  appropriate  soil  properties  and  pile 
boundary  conditions  were  provided  as  input  to  the  program.  The  results 
obtained  were  in  the  form  of  a  tabulation  of  horizontal  displacement, 
bending  moment,  sod  pressure,  and  modulus  of  horizontal  subgrade  reaction 
at  several  points  along  the  embedded  length  of  the  pile. 

Computer  Solutions.  The  data  measured  during  the  lateral  load  tests 
on  piles  consisted  of  lateral  displacements  and  slopes  of  the  pile  at  the  level 
of  the  ground  surface  for  various  magnitudes  of  later af  loads  Bending 
moments  as  a  result  ol  the  load  applications  above  the  ground  surface  were 
accounted  lor  in  the  computations  Computed  magnitudes  ol  displacement 
and  slope  am  presented  in  Futures  M  and  15.  respectively,  for  the  tests  at  El 
Centro  For  these  r omputahons.  v  lines  o*  k,  equa'  to  400S  and  200S  were 
ised  lor  the  do,  ‘o  ali<>n  and  the  Hooded  location,  respectively.  It  can  be 
sepn  that  the  (ompuled  results  are,  vmHi  few  exceptions,  somewhat  conserva 
Pec. 

Berause  the  previously  reported  results  of  tests  in  the  bay  mud  and 
the  hydruulii  hit  have  tx  en  referred  to  frequently  in  this  report,  a 
.  ornpansnn  o*  .jmp  uteri  and  measured  hun/or  tjl  displacements  during  tests 
•n  those  suns  is  provided  >n  f  igures  19  and  20.  respectively  Again,  the 
<  omputei  sol  it  ions  have  provded  sufln  ipntly  conservative  results 
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Horizontal  Load  on  Pile,  Q  (lb  x  10  ) 


Horizontal  Displacement  at  Groundline,  y  (in.) 


Horizontal  Displacement  at  Groundline, 


(in.) 


it  curves  from  pile  tests  in  bay  mud. 
represented  by  solid  lines.) 


Horizontal  Displacement  at  Groundline.  yg  (in.) 


Figure  20.  Results  of  lateral  load  tests  on  piles  in  the  hydraulic  fill. 


TABULAR  SOLUTION  FOR  LATERALLY  LOADED  PILES 

There  are  many  cases  where  it  is  desirable  to  have  a  reasonably 
accurate  solution  of  the  laterally  loaded  pile  problem  with  nonlinear  soil 
behavior  being  considered.,  but  in  most  of  those  cases  it  is  not  possible 
or  economical  to  obtain  access  to  a  computer  and  a  computer  program  of 
the  type  listed  in  the  Appendix  In  such  cases,  it  is  possible  to  utilize  a  more 
approximate  method  of  computation  which  has  been  developed  during  this 
study  This  method  utilizes  the  nonlinear  soil  pressure-displacement  relation¬ 
ships  presented  earlier,  while  computing  an  approximate  variation  of  the 
modulus  of  horizontal  subgrade  reaction,  k,  with  depth  These  magnitudes  of 
k  can  then  be  substituted  into  existing  design  charts,  and  pile  response  can 
be  estimated. 


35 


During  the  computer  study  discussed  earlier,  it  was  noted  that  in  all 
cases,  the  computed  variations  of  k  with  depth  could  he  approximated 
reasonably  with  the  weii-knowi.  function 

k  =  nh  Z  (10) 

where  nh  is  a  constant  of  horizontal  subyrade  reaction.  Equation  10  was 
reasonable  even  for  the  deposits  of  cohesive  soil  with  a  desiccated  surface 
layer.  Thus,  it  was  possible  to  develop  a  procedure  for  estimating  an  appro¬ 
priate  value  of  nh  for  a  given  pile  in  a  deposit  of  either  cohesive  or 
noncohesive  soil  for  a  given  condition  and  magnitude  of  loading.  It  is  noted 
tha*,  even  though  k  was  assumed  to  be  directly  proportional  to  depth, 
similar  design  procedures  can  easily  be  developed  with  other  assumed 
variations  such  as  step  functions  or  exponential  functions 

There  are  several  steps  involved  in  the  hand  solution,  and  these  steps 
are  described  before  they  are  presented  in  tabular  form. 

1.  Assume  a  trial  value  of  nh  for  the  soil.  It  has  been  found  that  by 
assuming  a  value  of  nh  larger  than  the  anticipated  actual  value,  convergence 
is  obtained  more  rapidly  than  by  assuming  lower  values.  Therefore,  an 
assumed  value  of  50  lb/in.3  tor  nh  is  recommended. 

2  Compute  a  trial  value  of  a  relative  stiffness  factor,  T,  defined  as 

follows 


w 

where  El  is  the  flexural  stiffness  of  the  pile 

3.  Determine  the  necessary  soil  properties  at  several  depths  along  the 
length  of  the  pile,  and  give  added  weight  to  the  soil  at  shallow  depths  in  a 
manner  similar  to  that  shown  later  in  Table  3.  For  cohesive  soils,  a  knowl¬ 
edge  of  the  undrained  shear  strength  at  those  depths  is  required.  For  granular 
soils,  the  angle  of  interna!  friction,  <p,  or  an  index  thereof  (Equation  7),  and 
the  effective  overburden  pressure,  pQ,  are  required. 

4.  With  the  values  of  El  and  nh  assumed  for  the  first  trial,  compute 
trial  values  of  pile  displacements  at  the  depths  for  which  soil  properties  were 
determined.  These  computations  can  be  made  with  the  use  of  procedures 
which  have  been  proposed  for  granu'ar  soils  8  The  resulting  design  charts 
have  been  presented  in  a  Navy  Design  Manual.1 7 
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5  With  the  trial  displacements  in  conjunction  with  the  measured  soil 
properties,  compute  a  corresponding  value  of  the  modulus  of  horizontal 
subgrade  reaction,  k,  at  each  depth  with  Equation  3.  For  cohesive  soils, 
Equations  5  and  6  can  be  used;  and  for  granular  soils.  Equations  7,  8,  and  9 
are  applicable. 

6.  A  plot  of  k  versus  depth  can  be  made  to  determine  the  next  trial 
value  of  nh  or,  if  necessary,  another  more  appropriate  distribution  of  k  with 
depth  can  be  determined.  In  the  tabular  solution,  a  procedure  is  presented 
for  computing  the  next  trial  value  of  nh  without  the  necessity  of  making  a 
plot  of  k  versus  2. 

7.  Using  the  next  trial  value  of  nh,  repeat  steps  2  through  6.  Perform 
a  sufficient  number  of  iterations  that  the  resulting  displacements  are  approx¬ 
imately  the  same  from  one  trial  to  another. 

8.  With  the  final  trial  value  of  nh,  compute  such  other  pile  response 
parameters  as  are  necessary  by  using  the  available  design  charts.1 7 

For  the  tabular  solution,  an  expediency  was  developed  to  allow  the 
computation  of  trial  average  values  of  nh  (step  no.  6)  without  the  require¬ 
ment  of  a  plot  of  k  versus  Z.  It  was  recognized  that  the  computed  values  of  k 
for  the  zone  of  soil  near  the  ground  surface  would  have  more  influence  on 
the  behavior  of  a  laterally  loaded  pile  than  would  those  values  for  depths 
near  the  point  of  inflection  of  the  displaced  pile.  It  was  noted  also  that  the 
point  of  inflection  typically  occurs  at  a  depth  equivalent  to  slightly  less  than 
21.  Therefore,  depths  for  which  computations  of  k  are  to  be  made  were 
defined  in  terms  of  T,  and  the  concentration  of  depths  for  computations  was 
specified  to  diminish  with  depth.  At  each  depth  a  corresponding  value  of  nh 
can  be  computed  with  the  trial  value  of  k  and  with  Equation  10;  a  simple 
arithmetic  mean  of  the  resulting  values  of  nh  provides  the  desired  weighted 
average 

Initially,  reasonable  values  of  depth  in  terms  of  T  were  determined 
intuitively.  Trial  solutions  then  were  compared  with  the  more  exact 
computer  solutions  using  the  program  presented  in  the  Appendix.  Minor 
adjustments  to  the  specified  values  of  depth  were  made  on  the  basis  of  the 
resulting  comparisons,  and  the  established  values  are  presented  in  Table  3. 
The  values  of  Ay  and  By  were  taken  from  design  curves  in  the  Navy  Design 
Manual1 7  for  an  infinitely  long  pile.  Therefore,  if  those  values  of  Ay  and  By 
are  used,  each  trial  solution  should  be  checked  to  ensure  that  the  pile  length 
is  greater  than  4T.8  It  was  found  from  this  study  that  the  proposed  tabular 
solution  resulted  in  computed  pile  response  data  with  sufficient  accuracy,  in 
comparison  with  the  computer  solutions,  that  the  computer  solutions  can  be 
supplanted  in  most  typical  cases. 
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Table  3.  Tabular  Form  of  Computational  Procedure 
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It  is  noted  that  each  tabular  solution  provides  a  value  of  nh  which  is 
applicable  only  for  the  specific  loading  condition  considered.  For  other 
loading  conditions,  a  different  value  of  nh  would  result  if  all  other  pile  and 
soil  conditions  remained  the  same.  This  behavior  (as  a  result  of  the  nonlinear 
soil  response)  is  demonstrated  in  the  following  example  problem. 

Example  Problem 

In  this  example  problem,  an  actual  field  test  is  utilized  in  order  that  a 
comparison  between  computer  results,  results  of  tabular  solutions,  and 
measured  test  data  are  available.  Test  No.  P4  is  considered.  The  pile  proper¬ 
ties  are  as  follows- 

Width,  B  =  16  in. 

Stiffness,  El  =  1 .69  x  101 0  lb-in.2 

Length,  L  =  324  in. 

A  profile  of  vane  resistance  versus  depth  was  presented  in  Figure  10  and, 
because  the  test  was  performed  in  the  dry  location  at  El  Centro,  the  larger 
values  of  S  are  appropriate  Two  levels  of  horizontal  load,  10,000  and 
15,000  pounds  are  considered,  along  with  the  moments  resulting  from  the 
31 -inch-height  of  the  point  of  application  of  the  loads. 

Because  a  comparison  of  theoretical  and  experimental  results  is  being 
made,  values  of  k|  =  4Q0S  are  used  in  both  the  tabular  solution  and  the 
computer  solution  instead  of  the  suggested  design  value  of  100S  indicated  in 
Equation  5.  In  both  cases,  the  values  of  p,  are  predicted  with  Equation  6. 

It  was  mentioned  during  the  preceding  discussion  of  the  tabular 
solution  that  relatively  rapid  convergence  can  be  achieved  by  assuming  an 
initial  trial  value  of  nh  greater  than  the  anticipated  actual  value.  It  has  been 
found  that  essentially  the  same  final  solutions  result  regardless  of  the  initial 
assumption  of  nh  ;  therefore,  a  low  value  of  nh ,  10  lb/in.3,  is  used  in  this 
example  problem  in  order  to  provide  a  demonstration  of  the  relatively  rapid 
convergence  resulting  even  with  a  poor  assumption  of  nh. 

The  entire  tabular  solution  is  presented  in  Table  4.  Nine  iterations 
were  performed  although  there  was  very  little  change  in  the  computed  dis¬ 
placements  after  the  fifth  iteration  The  resulting  value  of  0.18  inch  for  the 
displacement  at  the  groundline  is  compared  in  T able  5  with  the  corre¬ 
sponding  value  computed  with  the  computer  program  listed  in  the 
Appendix.  The  displacement  measured  during  the  field  test  is  presented  also.  It 
is  pointed  out  that  the  comparison  between  the  two  computed  values  is  not  as 
good  as  is  usually  found;  in  that  regard,  the  percentage  difference  between  the 
two  computed  values  for  a  load  level  of  1 5,000  pounds  is  somewhat  better . 
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Table  4.  Tabular  Form  of  Computational  Procedure  (Contd) 
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Table  4.  Tabular  Form  of  Computational  Procedure  (Contd) 
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Table  4.  Tabular  Form  of  Computational  Procedure  (Contd) 
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Table  5.  Comparison  of  Computed  Displacements 


i 

! 


Displacement  at  Groundline 
(in.) 

1 

|  Tabular 

Computer 

Test 

Solution 

Solution 

Data 

| 

0.18 

0.14 

0.15 

0.36 

0.30 

0.28 

Load 

(lb) 


10,000 

15,000 


The  realistic  handling  of  the  nonlinear  soii  response  is  evidenced  by 
the  fact  that  the  displacements  computed  by  the  tabular  method  are 
increased  100%  by  only  a  50%  increase  in  the  magnitude  of  load  from 
10,000  to  15,000  pounds.  Further  evidence  is  given  by  the  resulting  magni¬ 
tudes  of  nh  for  the  two  load  levels.  From  Table  4  it  can  be  seen  that  the 
computed  average  values  of  nh  for  a  10,000-pound  load  was  approximately 
91  lb/in.3  A  separate  solution  (not  presented)  for  a  load  of  15,000  pounds 
yielded  a  value  of  approximately  58  lb/in.3.  Therefore,  it  is  shown  that  the 
constant  of  horizontal  subgrade  reaction,  nh,  of  a  given  soil  deposit  is  not 
actually  a  constant. 

It  is  noted  that  if  estimates  of  other  pile  response  parameters  such  as 
slopes,  shears,  and  bending  moments  are  desired,  these  parameters  can  be 
obtained  by  using  the  available  design  charts1 7  and  the  values  of  nh  obtained 
during  the  displacement  computations. 


CONCLUSIONS 

From  the  results  of  lateral  plate  bearing  tests,  lateral  load  tests  on 
rigid  poles,  segmental  pile  tests,  lateral  load  tests  on  flexible  piles,  corre¬ 
sponding  soil  investigations,  and  theoretical  studies  of  lateral  soil  pressure- 
displacement  relationships,  the  following  conclusions  were  reached: 

1 .  For  the  design  of  laterally  loaded,  soil-supported  structures,  it  is  necessary 
to  account  for  the  nonlinear  response  of  the  soil  under  loading. 
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2.  Lateral  soil  pressure-displacement  cu'ves  can  be  determined  with 
sufficient  ease  and  accuracy  by  a  representation  of  these  relationships  with 
rectangular  hyperbolas  of  the  form 


pB  kj  +  pf  (  B  / 

3.  The  terms  kj  and  pf  in  the  above  relationship  can  be  defined  with 
sufficient  accuracy  in  terms  of  the  shear  strength  Cundrained  for  cohesive 
soils),  the  depth  below  the  ground  surface,  and  the  pile  width  for  both 
cohesive  and  noncohesive  soils.  With  noncohesive  soils,  a  knowledge  of  the 
effective  overburden  stress  versus  depth  is  required  also.  Expressions  for  kj 
and  p,  are  presented. 

4  These  lateral  soil  pressure-displacement  relationships  can  be  utilized  in  the 
computer  program  presented  in  the  Appendix  to  compute  the  response  of 
vertical  piles  to  lateral  loading.  The  computer  program  provides  a  finite 
difference  solution  in  an  iterative  form  in  order  that  the  nonlinear  soil 
behavior  can  be  considered. 

5.  For  cases  where  an  electronic  computer  is  not  available,  the  tabular  solu¬ 
tion  presented  is  an  expedient  means  of  obtaining  solutions  to  the  laterally 
loaded  pile  problem  while  considering  the  nonlinear  soil  response. 

6.  The  proposed  procedures  are  applicable  for  conditions  of  sustained 
loading.  For  cyclic  or  repetitive  loadings,  it  is  concluded  that  displacements 
and  slopes  should  be  increased  by  50%  over  the  computed  values  and  that 
bending  moments,  shears,  and  computed  soil  pressures  should  remain 
unchanged. 
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Appendix 


COMPUTER  PROGRAM 


A  format  for  a  once-through  solution  of  Equation  1  has  been 
suggested.8  This  format  was  followed  during  this  study,  but  it  was  not 
possible  to  perform  a  closed,  once-through  solution  because  of  the  nonlinear 
soil  response.  Therefore,  an  iterative  procedure  was  utilized. 

The  computer  code  for  an  iterative  finite  difference  solution  of 
Equation  2  is  presented  in  Table  A-t.  A  similar  computer  code  was  presented 
and  discussed  in  Reference  4,  but  the  code  has  been  changed  slightly  in  order 
to  provide  compatibility  with  the  more  appropriate  procedures  for 
predicting  soil  response.  In  addition,  other  minor  refinements  have  been 
made,  including  provisions  for  considering  tapered  piles  and  layered  soil 
systems. 

Comment  statements  are  included  in  the  computer  code;  however,  it 
is  desirable  to  have  some  additional  explanation  in  order  to  avoid  confusion, 
especially  with  regard  to  the  input  data.  All  of  the  input  data  are  arranged  in 
the  same  format  and  each  input  data  card  is  discussed  separately. 

Card  No.  1 

B  =  pile  width  at  the  groundline  (in.) 

El  =  flexural  stiffness  of  the 

pile  at  the  groundline  (lb-in.2) 

DX  =  incremental  lengths  of  the  pile 
or  pile  length  divided  by  53  (in.) 

PG  =  groundline  shear,  Qg  (lb) 

BMG  =  bending  moment  at  the  groundline,  Mg  (in. -lb) 

Card  No.  2 

CEI  =  variation  of  El  with  depth; 

EIA  (value  of  El  at  a  depth  Z) 

=  (CEI)  Z+  El  (lb-in.) 

T  =  number  of  iterations  required 

(1 1  iterations  were  sufficient  in  all  cases) 

BS  =  variation  of  pile  width,  B,  with  depth 
for  a  tapered  pile;  BA  (value  of  B  at  a 
depth  Z)  =  (BS)  Z  +  B  (in./in.) 
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Card  No.  3 


C(1)  through  C(5)  =  layer  designations  with  the  larger 
numbers  denoting  deeper  layers;  the  designations  consist  of 
the  numeral  1  for  clay  and  2  for  sand 

Card  No.  4 

ZC(1)  through  ZC(5)  =  depth  to  the  top  of  each  layer; 

ZC{1)  =  0  (in.) 

Card  No.  5 

This  card  contains  coefficients  for  the  computation  of  the  initial 
subgrade  modulus,  kj,  and  the  lateral  bearing  capacity,  pf, 
of  the  soil.  These  coefficients  were  varied  during  the  study, 
but  a  firm  set  of  values  has  now  been  proposed  as  follows: 

PFB  =  0.25 
CKIB=  100 
SKIA  =  2,000  (psi) 

AX  =  0.5 

Card  No.  6 

POI  =  effective  overburden  stress  at  the  ground  surface  (psi) 

P02  =  effective  overburden  stress  at  some  depth  Z2  (psi) 

Z2  =  depth  of  a  distinct  break  in  the  curve  of  pD  versus 
Z;  usually  the  depth  of  the  water  table  (in.) 

SL2  =  variation  of  p0  with  depth  below  the  depth  Z2; 

PO  (at  any  depth  Z  greater  than  Z2)  =  P02  +  SL2(Z  -  Z2) 

(psi/in.) 

Each  of  these  parameters  should  be  nonzero  in  order  to  avoid  having  zero 
denominators;  small  decimals  can  be  used  in  place  of  zeros. 
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Cards  No.  7  and  8 

SM{  1 )  through  SM(5)  =  values  of  vane  resistance,  S,  at  the 
corresponding  depths  of  ZM(1)  through  ZM(5);  again,  the 
values  should  be  nonzero  !psi) 

SMB  =  variation  of  vane  resistance  with  depth  below  a  depth  of 
ZM{5);  S  [at  any  depth  Z  greater  than  ZM(5)] 

=  SMB  [Z  -  ZM(5)]  +  SM{5)  (psi/in.) 

ZM{  1 )  through  ZM(5)  =  corresponding  values  of  depth, 

ZM{1)  =  0 

A  brief  explanation  of  the  procedure  used  for  entering  vane  resis¬ 
tances  into  the  computer  solution  is  in  order.  As  was  mentioned  previously, 
the  behavior  of  laterally  loaded  piles  is  influenced  most  extensively  by  the 
properties  of  the  soil  at  shallow  depths.  In  addition,  the  vane  resistance  of 
soil  at  depths  below  the  surface  layers  often  increases  approximately 
uniformly  with  depth.  Therefore,  the  typical  procedure  used  in  representing 
S  versus  Z  was  to  use  SM(1 )  through  SM(5)  to  provide  a  fairly  accurate 
representation  at  shallow  depths  and  to  use  SMB  for  a  less  exact  average 
representation  at  greater  depths. 

The  computer  code  provides  a  trial  solution  for  initiation  of  the 
iterative  solution  by  using  an  arbitrary  value  of  k  of  300  psi  at  all  points 
along  the  length  of  the  pile.  With  the  resulting  trial  value  of  displacement  at 
each  point  and  the  soil  conditions  provided,  a  new  value  of  k  is  computed 
for  each  point.  With  the  revised  values  of  k,  the  finite  difference  solution  is 
repeated  and  new  displacements  are  obtained.  This  procedure  is  repeated 
until  convergence  is  achieved.  It  was  found  that  satisfactory  convergence 
occurred  in  all  cases  with  1 1  iterations,  and  the  entire  solution  required  less 
than  2  minutes  of  time  on  an  IBM  1620  computer 

The  accuracy  of  the  computer  program  was  checked  in  two  ways. 
First,  bending  moments  computed  at  the  ground  elevation  were  compared 
with  magnitudes  of  bending  moment  used  as  input  to  the  solutions  There 
were  negligible  differences  in  the  magnitudes,  so  it  was  concluded  that  the 
counters  and  index  parameters  were  computed  properly.  The  second  method 
of  checking  the  program  consisted  of  obtaining  computer  solutions,  with  k 
assumed  constant  with  depth  and  displacement,  and  then  comparing  these 
solutions  with  results  of  a  theoretical  study  using  an  analog  computer  Again, 
the  comparison  was  favoiable  It  was  concluded  that  the  difference-equation 
method  of  solution  with  53  increments  considered  along  the  pile  length 
provided  results  of  sufficient  accuracy  when  realistic  soil  parameters  were 
considered. 
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Table  A  l.  Computer  Code 


DIMENSION  Htlll).  G157).  Y(59),  SMt6>»  2M(6).  Ct6),  ZC(6) 

DATA  INPUT 

PILE  WIDTH»  PILE  STIFFNESS.  EFFECTIVE  LENGTH  DIV.  8Y  53.  GROUNDLINE 
SHEAR.  GROUNDLINE  MOMENT. 

5  READ  100.8.EI.DX.PG.3MG 

STIFFNESS  VARIATION.  ITERATIONS.  WIDTH  VARIATION. 

READ  100.CEI.T.8S 
C  SAND  OR  CLAY  LAYER. 

READ  100<C(1)<C(2),C(3),C(4),C(5> 

C  DEPTH  TO  TOP  OF  LAYER. 

READ  100.  ZC( 1 ) »  2CI21.  2C(3).  ZC(4),  2C(5) 

C  SOIL  PRESSURE  -  DEFLECTION  CONSTANTS. 

READ  l^O,  PF8.CKI8.SK  T  A . AX 

C  EFFECTIVE  OVERBURDEN  PRESSURE  AND  WATFR  CONDITIONS  FOR  SAND. 

43  READ  100.P01.P02.Z2.SL2 

C  UNDRAINED  SHEAR  STRENGTH  OATA  t  VANE  SHEAR  I. 

33  READ  lOO.SMt 1) ,SM<2) ,SM(3) ,SM(4> «SM(5> .SM8 

C  OEPTHS  FOR  SHEAR  STRENGTH  OATA. 

READ  lOO.ZMtl) .ZM(2) ,2M(3> .ZM<4> .ZMI5) 

100  FORMAT (5E15.8.F5.2) 

PUNCH  100,3. El. DX.PG.8HG 

PUNCH  100.CEI.T.3S 

PUNCH  100,CUl.C(2t.C(3).C(4j,Ct5> 

PUNCH  100.  2CHJ.  ZC  <  2 !  •  ZC(3).  ZC(4),  ZCI5) 

PUNCH  100,  PFB.CKIB.SKIA.AX 
PUNCH  100.  P01.  P02.  22.  SL2 
PUNCH  100,SK(1).SM121 «SM< 3 > »SM<4) ,SM( 51 .SMB 
PUNCH  lOO.ZMtl), ZM<2J ,ZM< 3 > .2Mt4 ) ,ZM[ 5 > 

C  ITERATIONS. 

M=T 

C  DIFFERENTIAL  increments. 

N=S3 

DX2=DX*DX 

DO  55  L=1.M 

SET=L-1 

N=N.  1 

IJ  =  5 

IK-5 

IF(SET)24.34,24 

34  DO  44  J=1 ,N 
I  =J+3 

C  INITIATION  OF  SOIL  PARAMETERS. 

44  G(I)=300.»DX2*DX2/EI 
GO  TO  45 

C  START  AT  TIP  OF  PILE. 

24  DO  56  J«1,N 
I  =  J+3 
EYE=N-J 
Z=EYE»DX 

C  ACTUAL  PILE  WIDTH. 

BA=eS*Z*B 

C  ACTUAL  PILE  STIFFNESS. 

EIA=CEI*Z+EI 

C  SHEAR  STRENGTH  CALCULATIONS. 

I F ( Z-ZM( 51)76,86 ,86 

86  S=SMB*(Z-ZM(5))+SM(5> 

GO  TO  83 

76  IF(Z“ZM( IJ-1 ) )79,89,8R 
79  I J= I J-l 
GO  TO  76 

89  S=(SM( I Jl-SMt I J-l ) >*<Z-ZM( IJ-1 ) ) / ( ZM ( I J ) -ZM ( I J-l ) )+SM( I J-l ) 

83  IF(Z-ZC< IK))87, 97,97 

87  IK= I K-l 
GO  TO  83 

C  IS  LAYER  SAND  OR  CLAY. 


Continued 
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T able  A- 1 .  Computer  Code  (Contci) 


97  K !ND=C ( IK ) 

GO  TO  1 69,93 ) .KIND 

C  HYPERBOLIC  PARAMETERS  FOR  CLAY. 

69  PF  =  PFb+S»Z/<8*»AX) 

CKI  =  CK  IB  «  S 
GO  TO  65 

C  HYPERBOLIC  PARAMETERS  FOR  SAND. 

93  IF ( 2-221 37. 37,47 
37  PO=P01+(P02-P01 )*2/Z2 
GO  TO  35 

47  P0=P02+SL2*<Z-Z2l 

35  PHI  =  ( 25.4-PO/ <2.94+0. 1225+PO) I *S**0. 1 6 

HNO=0.664»< 10.»*(0.038»PH! I 1+0.049* t 10.»* < 0.034+PHI ))*Z/BA 

PF=HNQ#PO 

CKI  ■  SKIA 

C  RECTANGUL AR-HYPER80LA  RELATIONSHIP. 

65  YOP=6A/CKl+Y< I )/PF 
CNG=DX2*DX2»BA/£ I A 
56  G( I )=CNG/YOP 

C  DIFFERENCE  -  EOUATION  SOLUTION. 

45  H.'4)»2./'2.+G<4)  ) 

H(  5 ) =2.*H( 4 ) 

H(6!=l./<5.+G<5)-2.»H<5)) 

H(7)=H(6)*(4.-H(5) ) 

N=N-2 

DO  JO  1=1, N 
J*2»( I+5I-4 
K  =  1+5 

H(J)=l./(6.+6(K|-H( J-41-HI J-l I*f4.-H( J-3II ) 

10  H(J+1)=H(J1»(4.-H(J-1  II 
CJ4=2.*DX2*DX*PG/EI 
CJ3=0X2*BMG/E I 
N=N+1 
K=2»N+4 

Y(N)=(CJ4*H(K)*( 1 «-H ( K-2 I ) I 

Y(N)-=Y(N|-CJ3«H(K)*(2.  +  2.*H(K-2»-HtK-2)*H(K-3M 
Y(N|=Y(N)+CJ3*H(K+1 I 

A=4.-4.*H(K-l)+4.»H(K-2)-H(K-4>+H(K-l >*H(K-3! 

A=A-2.*H ( K-2 )+H(K-3 I +H( K-2 1 »H<  K-4) 

Y(N+4)=Y<N)/< l.-H(K-2)-2.*H<K+l >+H<K-l )*H(K+1 )+HCK)*A) 
Y(N+5)=(CJ3+Y(N+4)*(2.-H(K-l) ) I /(l. -HI K-2) ) 

Y(N+6)  =  (H(K+l)*Y(N+5)-Y(fl+4)  )/H(KI 

DO  20  J= 1 *N 

I=N-J+4 

K=2*I-4 

20  Y( f)=-H(K)*Y( I+2)+H(K+l )*Y( 1+1 ) 

Y(3)=-Y(5)+2.*Y(4) 

Y(2)=Y(6)-2.*Y(51+2.»Y(31 
SL= ( -Y  <N+5 1 +Y  < N+3 11/(2. *0X1 
55  PUNCH  100* Y ( N+4  J  *  SL 
N=N  +  1 

DO  80  J  =  1 «N  *  2 
I =N— J+4 
EYE=J-1 
Z  =  EYE  «  DX 
BA  =  BS  »  Z  +B 
EIA  =  CF!  *  Z  +FI 

BM=EIA*( Y( 1+1 )-2.*Y( I)+Y( 1-1 ) ) /DX2 

CNG=DX2»DX2»BA/EIA 

P=Gt I i*Y< I ) /CNG 

CK=P*BA/Y ( I ) 

C  OEPTH,  DEFLECTION,  MOMENT,  SOIL  PRESSURE,  SUBGRADE  MODULUS. 

80  PUNCH  100,  Z,  Y(I),  BM,  P,  CK 
GO  TO  5 
70  PAUSE 
END 
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LIST  OF  SYMBOLS 


Ay  Parameter  from  design  curves 

B  Width  of  pile,  in. 

By  Parameter  from  design  curves 

d  Horizontal  displacement  of  test  segment,  in. 

El  Flexural  stiffness  of  pile.  lb-m  2 

H  Height  of  test  segment,  in. 

h  Height  above  ground  surface  of  lateral  load 
on  pile,  in. 

k  Modulus  of  horizontal  subgrade  reaction,  psi 

kj  Value  of  k  at  zero  displacement,  psi 

k0  Value  of  k  assumed  constant  with  depth,  psi 

L  Embedded  length  of  piie,  in. 

Lj  Liquidity  index  of  soil,  % 

M  Bending  moment  in  pile,  in.-Sb 

Mg  Applied  bending  moment  at  ground  level,  in. -lb 

N  Standard  penetration  resistance 

of  soil,  blows/ft 

n^  Constant  of  horizontal  subgrade 

reaction,  lb/in.3 

p  Horizontal  soil  pressure  psi 

Pf  Horizontal  bearing  capacity  of  soil,  psi 


p0  Vertical  effective  stress  due 

to  overburden,  psi 

Q  Horizontal  load  on  pile,  lb 

Qg  Magnitude  of  y  at  ground  level,  in. 

R  Relative  stiffness  factor  for 

k  constant  with  depth,  in. 

S  Vane  resistance  of  soil,  psi 

T  Relative  stiffness  factor 

for  k  directly  proportional  to  depth 

y  Horizontal  displacement  of  pile.  in. 

yg  Magnitude  of  y  at  ground  level,  in. 

ym  Magnitude  of  y  as  a  result  of  a 

bending  moment  loading 

Yq  Magnitude  of  y  as  a 

result  of  a  shear  loading 

Z  Depth  below  ground  level,  in. 

a  Shape  factor 

Og  Slope  of  pile  at  ground  level 

0  Angle  of  internal  friction 

of  the  sand,  degrees 
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1.  Please  make  the  following  pen  and  ink  corrections: 
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absolute  value  of  displacement,  lyl,  should  be  substituted  for  y. 


PETER  D.  TRIEM 
By  direction 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

Springfield,  V«.  22151 


1 1  NO  .NCEL. 10440/1  (REV.  4.40) 


